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Abstract 

Wireless power transfer (WPT) prolongs the lifetime of wireless sensor network by providing 
sustainable power supply to the distributed sensor nodes (SNs) via electromagnetic waves. To improve 
the energy transfer efficiency in a large WPT system, this paper proposes an adaptively directional 
WPT (AD-WPT) scheme, where the power beacons (PBs) adapt the energy beamforming strategy to 
SNs’ locations by concentrating the transmit power on the nearby SNs within the efficient charging 
radius. With the aid of stochastic geometry, we derive the closed-form expressions of the distribution 
metrics of the aggregate received power at a typical SN and further approximate the complementary 
cumulative distribution function using Gamma distribution with second-order moment matching. To 
design the charging radius for the optimal AD-WPT operation, we exploit the tradeoff between the 
power intensity of the energy beams and the number of SNs to be charged. Depending on different 
SN task requirements, the optimal AD-WPT can maximize the average received power or the active 
probability of the SNs, respectively. It is shown that both the maximized average received power and the 
maximized sensor active probability increase with the increased deployment density and transmit power 
of the PBs, and decrease with the increased density of the SNs and the energy beamwidth. Finally, we 
show that the optimal AD-WPT can significantly improve the energy transfer efficiency compared with 
the traditional omnidirectional WPT. 
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I. Introduction 

Wireless sensor networks (WSNs) eonsist of small-size, low-power and distributed sensor 
nodes (SNs) to monitor physical or environmental conditions [HI. WSNs are often required to 
operate for long periods of time, but the network lifetime is constrained by the limited battery 
capacity and costly battery replacement at SNs. To extend the network lifetime, it is desirable 
to recharge the SNs in an undisruptive and energy efficient way. 

RF-enabled wireless power transfer (WPT) [|2l provides a controllable and sustainable power 
supply to sensor network by charging SNs via electromagnetic (EM) waves O-llSl. There are 
mainly two types of WPT: omnidirectional WPT and directional WPT. For omnidirectional 
WPT, the energy transmitter or so-called power beacon (PB) broadcasts the EM waves equally 
in all directions regardless of the locations of the energy receivers. According to the law of 
conservation of energy, the energy radiated in the direction of energy receivers accounts for only 
a small fraction of the total radiated power. Since the EM waves fade rapidly over distance, it 
may require excessively high transmit power to charge an energy receiver via omnidirectional 
WPT, which may not be energy efficient. In contrast, for directional WPT with antenna arrays, 
the PB concentrates the radiated energy in the directions of the energy receivers, i.e., via energy 
beamforming, which enhances the power intensity in the intended directions. The energy transfer 
efficiency is thus improved with the consequent reduction of transmit power to reach the target 
received power. 

Most of the literature on directional WPT (see [|2]| and references therein) has focused on point- 
to-point and point-to-multipoint transmissions. For a large-scale WSN, the SNs are often in large 
quantities and are usually distributed with random locations. There are two main challenges in 
the design of directional WPT for a large-scale network. On the PB-side, it is challenging to 
adapt the energy beamforming strategy to the random locations of the SNs, e.g., to decide which 
SNs to serve, how many beams to generate and the beamwidth of each beam, etc. On the SN- 
side, it is difficult to analyze the aggregate received power from a large number of PBs in the 
network, where the radiation directions and energy intensity may vary for each PB. 

In this paper, we aim at tackling the above two challenges. The paper structure and main 
contributions are given as follows. 

• Energy-efficient AD-WPT scheme to power a large-scale sensor network: To address the 
PB-side challenge, we propose an adaptively directional WPT (AD-WPT) scheme in a large- 
scale sensor network in Section II, where the energy beamforming strategy of the PBs is 
adaptive to the nearby SN locations that are within the energy-efficient charging radius. To 
deal with the tradeoff between the power intensity of the energy beams and the number 
of SNs served by each PB, we design the charging radius to achieve optimal AD-WPT 
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for different performance targets, i.e., average power maximization or active probability 
maximization. 

. Analysis of harvested power using stochastic geometry. To address the SN-side challenge, in 
Section III, we successfully derive the closed-form expressions of the distribution metrics, 
e.g., Laplace transform, mean and variance, of the aggregate received power at a typical 
SN from the large-scale PB network using the tools of stochastic geometry |l6l-|[8l. The 
complementary cumulative distribution function (CCDF) of the received power is also 
analyzed. As it is difficult to obtain the analytical CCDF expression, we further approximate 
it using Gamma distribution with second-order moment matching. 

• Optimal AD-WPT for average power maximization: In flexible-task WSN, the SNs operate 
in a cooperative manner on power adaptive sensing tasks. To achieve the optimal AD-WPT, 
we design the optimal charging radius to maximize the average received power of the 
SNs in Section IV. We show that the maximized average received power increases with the 
increased PB power and density, while it decreases with the increased energy beamwidth and 
SN density. In addition, the optimal AD-WPT greatly improves the average received power 
compared with the traditional omnidirectional WPT, especially when PB power/density is 
high. 

• Optimal AD-WPT for active probability maximization: In equal-task WSN, the SNs operate 
in an independent manner on equal quantity of sensing tasks, where an SN is active if its 
received power is larger than the operational power threshold. To achieve the optimal AD- 
WPT, we design the optimal charging radius to maximize the active probability of the SNs 
in Section V. It shows that the optimal AD-WPT can enhance the sensor active probability 
compared with omnidirectional WPT, especially when the PB power/density is not high. 

In Section VI, the numerical results are shown and discussed. Finally, conclusions are drawn in 
Section VII. 

A. Related Literature 

Omnidirectional WPT has been studied recently in [|9ll- [fT^ . In dU, a point-to-point omni¬ 
directional WPT is investigated, where the receiver utilizes part of the harvested energy for 
decoding the information in the received signal. In ifTOl . the downlink energy transfer in a 
broadcast network is studied for throughput maximization. In iHTI . a stochastic geometry based 
model is considered for a cognitive radio network, where the secondary transmitters harvest RF 
energy from the nearby primary transmitters. IfT^ investigates the downlink energy transfer in 
a large-scale wireless network by considering finite and infinite battery capacity. 

The directional WPT has been addressed in [[T3ll - [[T^ . In [[T3ll . energy beamforming is studied 
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Fig. 1. System model of AD-WPT (illustrative example of N = 4). The circular areas with radius p are the charging regions 
of the PBs. The shaded sectors in the charging regions are the active sectors of the PBs. 


in a broadcast network where the transmitter steers the energy beams towards the reeeivers to 
maximize their received power. In [fT4]| and [fTSll . energy beamforming is designed in a MIMO 
broadeast network jointly with information beamforming, where the transmitter adjusts the beam 
weights to maximize the reeeived power and information rate at different receivers. In ITT^ . each 
mobile node in a cellular network is charged by its nearest PB via energy beamforming. For 
the simplieity of analysis, only the reeeived power from the nearest PB is eonsidered and the 
received energy from all other PBs is omitted in ifT^ . 

To the best of our knowledge, this paper is the first study of direetional WPT by using adaptive 
energy beamforming for a large-seale network and the resulting aggregate received power from 
all PBs with AD-WPT is rigorously eharaeterized. With the proposed AD-WPT seheme, the 
energy transfer efficieney in the large-seale network can be greatly enhanced eompared with the 
traditional omnidireetional WPT. 


II. System Model 

We eonsider a wireless eharging network as shown in Fig. 1, where a PB network wirelessly 
eharges an SN network via energy beamforming. Eaeh PB radiates EM waves with wavelength 
u using transmit power Pp. The PBs and SNs follow two independent homogeneous Poisson 
Point Proeesses (PPPs) ^p = {Xi} and = {Yj} with density Xp and A^, respeetively, where 
Xi and Yj represent the eoordinates of the PBs and SNs in plane. 

In the following, we first propose a power transfer seheme with adaptive energy beamforming 
and then diseuss the power intensity in the direetions of the energy beams. 
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A. AD-WPT Scheme 

Due to the fast attenuation of the radio power over the distanee, it is more energy effieient 
for the PBs to foeus the energy to eharge the nearby SNs. With antenna arrays, a PB is able to 
form an energy beam in a eertain direetion or generate multiple beams simultaneously towards 
different direetions ifTTl . In this subseetion, we propose an AD-WPT seheme where the PBs 
adapt the beamforming strategy to the random loeations of the SNs. 

To deeide whieh SNs to eharge, we define charging region as a eireular region eentered at 
each PB with charging radius p, as shown in Fig. 1. Each charging region is divided into N 
equal sectors Ci, • • • , Cat, where N is usually a small positive integer due to physical constraint 
of antenna design. We consider that a PB is aware of the existence of the SNs inside each of its 
sectors, e.g., via the SN feedback over control channels. A sector is considered to be active if at 
least one SN falls into this sector. Denote M as the random number of active sectors of a PB, 
e.g., PBj, where 0 < M < N. The adaptive beamforming strategy of PBj is given as follows. 

• SN’s absence in charging region: If no sector of PBj is active (M = 0), PBj works as an 
omnidirectional antenna that radiates energy equally in all directions (to help power SNs 
outside the charging region). 

• SN’s presence in charging region: If at least one sector of PBj is active (M > 1), PBj 
generates M equal-power energy beams in the directions of the M active sectors. 

We use equal power allocation among the energy beams of a PB for the ease of analysis. In 
Section VI-C, we will show that equal power allocation is descent as compared with some other 
unequal allocation choices. 

From an SN’s point of view, the received power from the PBs is discussed as follows. 

• Inside charging region (or within radius p): An SN can be intentionally and efficiently 
charged by one or more PBs whose charging regions cover its location. 

. Outside charging region (or beyond radius p): When an SN is located outside the charging 
regions of the PBs, the SN still receives RF energy from the PBs if it is aligned with the 
energy radiation directions of the PBs. 

We further explain the proposed AD-WPT scheme with the example of = 4 in Fig. 1. It 
is observed that PBi detects three nearby sensors, i.e., SNi, SN2 and SN3, which fall into three 
out of four sectors of its charging region. As a result, PBi adaptively generates three energy 
beams in the directions of north-east, north-west and south-west to directionally charge the three 
sensors. At the same time, PB 2 detects three sensors, i.e., SN3, SN4 and SN5, which fall into 
two sectors of its charging region, and thus two adaptive energy beams are generated towards 
these SNs. In particular, notice that SN3, which is within the overlapping area of the charging 
regions of PBi and PB 2 , is thus intentionally charged by the two PBs at the same time. SNi is 
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intentionally charged by PBi while it also receives energy from PB 2 and PB 3 since its location 
is aligned with the south-west energy radiation directions the two PBs. 

B. Antenna Gain under AD-WPT 

When a PB is directional, the power intensity in the directions of energy beams improves 
compared with the case when the PB is omnidirectional. The ratio of power intensity between 
directional and omnidirectional antenna is defined as the gain of directional antenna G (G > 1 ) 
[fTTl . In the unintended directions of the directional PB, the power intensity is zero. In the 
following, we evaluate G given that M out of N sectors of the PB are active. 

If none of the sectors of a PB is active (M = 0), as discussed, the PB behaves as an 
omnidirectional antenna with the uniform gain in all directions, i.e., 

Gm = 1 , for M = 0 . ( 1 ) 

If M out of N sectors of the PB are active (M > 1), the PB forms M (M < N) energy beams 
with equal power in the direction of each beaml^rBy the law of conservation of energy, the total 
radiated power for directional and omnidirectional antenna is the same. Since the directional 
antenna concentrates the energy from the directions of N sectors into M sectors, the power 
intensity in the intended directions becomes N/M times of that of the omnidirectional antenna. 
Therefore, given M energy beams at the PB, the antenna gain in the direction of each energy 
beam is approximated as 

Gm = N/M, for M = !,■■■ ,N. (2) 

From (HI) and Q, we see that the proposed AD-WPT is equivalent to the omnidirectional WPT 
with uniform gain when M = 0 or M = N. 

The antenna gains and number of energy beams of the PBs are related to the charging radius 
p. As p —)■ 0, no SN is inside the charging regions (M = 0) and all PBs radiate energy in 
N directions with gain Gq = 1 as omnidirectional WPT. As the increase of p, more sectors 
of the PB are likely to be activated due to the increased number of SNs inside the charging 
region. The number of beams that most PB radiate with decreases from to 1 sharply and 
then increases from 1, 2, • • •, to A^. The corresponding antenna gain increases from Gq = 1 to 
Gi = A^ and then decreases from Gi = A^, G2 = y, ■ • •, to Gn = 1. As p — )■ cxd, AD-WPT 
is again equivalent to omnidirectional WPT with Gat = 1 in all A^ directions. As we can see, 
there is a tradeoff between the antenna gain Gm and the number of beams of the PBs. When 


'For simplicity, we assume the side lobes are negligible and the radiated energy is uniformly distributed across each energy 
beam. 
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the PB concentrates energy on fewer beams, the power intensity of each beam increases but at 
the cost of charging fewer SNs. To address the above tradeoff, the optimal charging radius is 
crucial in the AD-WPT design and will be analyzed in Section IV and Section V for different 
SN network requirements. 

III. Characterization of SNs’ Received Power Using Stochastic Geometry 

In this section, we first study the aggregate received power at a typical SN from all PBs and 
then use stochastic geometry to analyze the distribution of the received power. 

Consider a typical sensor node SNq at the origin and an arbitrary PBj at location Xj. If PBj 
radiates energy with gain Gm (for M = 0,1, • • • , iV) towards SNq, the received power at SNq 
from PBj is ifT^ 

P; = PpGMa [max (||X,||/do, 1)]”“ , (3) 

where Pp is the transmit power of PBj, a is the path loss exponent, a is a unitless constant 
depending on the receiver energy convention efficiency, antenna characteristics and average 
channel attenuation^ The Euclidian distance between PBj and SNq is represented by ||Xj||, 
and do is a reference distance for the antenna far field. The received power from each PB is 
taken by averaging over the short-term fading. We adopt the non-singular path loss model |I71 
to avoid [||Xj||/(io]~'^ > 1 for ||Xj|| < do. Without of the loss of generality, we use do = 1 
throughout the paper. 

Equation ([3]) holds if PBj radiates energy with gain Gm towards SNq, where Gm is given in 
([1]) or dH) depending on the number M of active sectors of PBj. By considering all PBs in the 
network, the aggregate received power at SNq is 

Ps = ^ P^l (SNq receives energy from PBj with Gm) ■ (4) 

The indicator function equals one if both the following conditions are satisfied: 

• Condition 1: PBj has M active sectors; 

• Condition 2: SNq is in one of the M radiation directions of PBj given PBj has M active 
sectors. 

We see that both conditions are related to the distance between SNq and PBj. If SNq is inside 
the charging region of PBj, PBj generates at least one beam towards SNq (M > 1). If SNq is 
outside the charging region of PBj, SNq may not be in the radiation direction of PBj and M 
may vary from 0 to X. 


^For empirical approximation, cr is sometimes set to free-space path loss at distance do assuming omnidirectional antennas, 
i.e., a = 201ogj^g dB flSl. where v is the wavelength of the radio waves. 
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According to the distance between PB* and SNq, we elassify the PBs into two groups: near 
PBs with ||Xj|| < p, and/ar PBs with ||Xj|| > p. We draw an equivalent eharging region 
eentered at SNq with radius p and denote 6(o, p) and 6(o, p) as the regions inside and outside 
this eharging region, respeetively. We define two indieator funetions 9^ and 9^ to deseribe the 
events that SNq reeeives power from the PB with Gm conditioned on this PB is a near PB or 
far PB, respeetively, i.e., 

9^ = 1 [SNq reeeives energy from PBj with Cm I ||-^i|| < p] (5) 

and 

9^=1 [SNq reeeives energy from PBj with Gm | ||-^i|| > p] , (6) 

where subseripts n and / denote the near and far PBs and superseript M denotes the number 
of aetive seetors of the PB. 

We denote Pg^n as the aggregate reeeived power from the near PBs and Pgj as the aggregate 
reeeived power from the far PBs that radiate energy towards SNq. By summing them up, we 
rewrite Pg as 

Ps = Ps,n + Psj^ (7) 

where 

Ps,n = PpCT GmC [max (||X,11,1)]-“ (8) 

Xi£^pf]b{o,p) 

and 

Psj = Pp^ G-M^f [max (||X,11,1 )]-“. (9) 

Xi£^pf]b{o,p) 

As a speeial ease of = 1, all PBs are omnidireetional radiators with gain of 1. The aggregate 
reeeived power at SNq from all omnidireetional PBs is 

pomm^p^^ ^ [max (||X,11,1)]-“. (10) 

To fully eharaeterize the reeeived power distribution, we usually use Laplaee transform, whieh 
is however, diffieult to be derived direetly from (|7]). In the eonditional events of 9^^ and 9^, 
the gain Gm of PBj is also related to the loeations of other nearby SNs of PBj whieh are 
unknown. Moreover, sinee Gm vary for eaeh PB, the PBs that radiate power with Gm towards 
SNq ean be regarded as a heterogeneous network for whieh the Laplaee transform is hard to 
eharaeterize. In the following diseussions, we use an alternative method by taking the privilege 
of the independent thinning |8l of the network. For the near PBs and the far PBs, respeetively, we 



8 

thin the heterogeneous network into multiple homogeneous networks with certain probabilities, 
where in each homogeneous network the PBs radiate energy towards SNq with the same gain 
Gm- We have M = 1, • • • , iV for the near PBs and M = 0,1, • • • , iV for the far PBs. After 
analyzing the Laplace transform of the received power distribution in each homogeneous network, 
we finally derive the distribution metrics of the aggregate received power from all PBs at SNq. 


A. Power Reception Probability given PB Location 

First, we derive the thinning probabilities of the near PBs and the far PBs. As discussed 
previously, SNq receives power from PBj with gain Gm if both Conditions 1 and 2 are satisfied. 
As for Condition 1, PBj transmits with gain Gm if it has M active sectors. We derive the active 
probability of each sector as follows. As SNs follow PPP with density As, the number of SNs 
inside a charging region is a Poisson random variable with mean AsTrp^. When the charging 
region is equally partitioned into N sectors, the number of SNs inside one of these N sectors 
is also a Poisson random variable, denoted by /, with mean /N, and the probability mass 
function is given by 

Pr {I = k) = exp (—XsTip'^/N) , k = 0,1, ■■ ■ (11) 

kI 

The probability that no SN is inside a sector is thus 

p = Pr (/ = 0) = exp (^—XsTip'^/N) . (12) 


Therefore, the active probability of a sector is the probability that at least one SN is inside this 
sector, which is given by 


q = l—p=l — exp [—\s7rp‘^/N) . 


(13) 


Denote and as the conditional probabilities that SNq receives energy from PB* with 
antenna gain Gm given PB* is a near PB and a far PB, respectively. Based on p, q. Conditions 
1 and 2, we derive p^ and p^ as follows. 

1) Near PBs: If ||Xj|| < p, PBj radiates energy in at least the direction towards SNq (M > 1). 
Condition 2 is thus satisfied. Given PBj is a near PB, the conditional probability that PB* radiates 
with gain Gm is 


N=(^ 1\ JV-M M-l 


(14) 


which is the probability that the rest M — 1 out of — 1 sectors of PB* have SNs. Given PBj 
is a near PB that radiates with gain Gm, the conditional probability that SNq receives energy 
from PBj is 




M 

n 


1 . 


(15) 
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Since rj^ = , we obtain the following lemma. 

Lemma 1: Given PBj is a near PB, the eonditional probability that SNq reeeives energy from 
PBj with gain Gm is 

= ( 16 ) 


2) Far PBs: If ||Xi|| > p, PB* may not radiate energy towards SNq (M = 0, • • • , iV). SNq 
reeeives energy PB* with Gm if both Conditions 1 and 2 are satisfied. Given PBj is a far PB, 
the eonditional probability that PBj radiates with gain Gm is 

(, for M = 0 (17a) 


u 


M 

f 




for M 


.N. 


(17b) 


Given PB* is a far PB that radiates with gain Gm, the eonditional probability that SNq reeeives 
energy from PB* is 



1 , 


for M 


M 

iV’ 


for M 


0 

I,--- ,iV. 


(18a) 

(18b) 


Sinee pY = we obtain the following lemma. 

Lemma 2: Given PBj is a far PB, the eonditional probability that SNq reeeives energy from 
PBj with gain is 



'N-1 

M-1 


pN-MqM^ 


for M 
for M 


0 

!,■■■ ,N. 


(19a) 

(19b) 


B. Characterization of Received Power via Laplace Transform 

In this subseetion, we derive the Laplaee transform of the distribution of Pg to eharaeterize 
the reeeived power at SNq. 

Define ^Y the set of PBs with gain Gm and as the set of PBs that radiate energy 
towards SNq. The set of near PBs within b{o, p) that radiate energy with gain Gm towards SNq 
is 

n n p). for Af=i, ■ ■ ■, v <20) 

whieh is obtained through the independent thinning |[8l of near PBs with new density XppY^ 
where pY given in Lemma 1. The near PBs ean be regarded as a heterogeneous network 
eonsisting of a group of homogeneous networks eaeh with antenna gain Gm and density XppY ■ 
Similarly, the set of far PBs within b(o, p) that radiate energy with gain Gm towards SNq is 


=< n n ^=0’ • • • ’ 


( 21 ) 
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which by the independent thinning of far PBs with new density Xprjf, where r]Y is given in 
Lemma 2. The far PBs that radiate power towards SNq can be regarded as another heterogeneous 
network consisting of a group of homogeneous networks each with gain Gm and density Xprjf. 
Note that SNq reeeives zero power from the far PBs that does not radiate energy towards SNq. 
In the following, we derive the Laplaee transform of the reeeived power distribution in eaeh 
homogeneous network, and then derive that of the aggregate reeeived power from all PBs. 

We rewrite the aggregate reeeived power at SNq from all the near PBs and far PBs in ([7]) as 

N N 

P. = P.,n + P.J = Gn + E (22) 

M=1 M = 0 


where 


P^n = Pp^ Gm [max (||X,11,1)]-“ (23) 

is the aggregate reeeived power from the near PBs with gain Gm and 


P^f = PpCT Y. Gm [max (||X, 11 , 1 )]-“ 




M 

PJ 


(24) 


is the aggregate reeeived power from the far PBs with gain Gm- Sinee we adopt the non-singular 
path loss function [max (|| 2 fj||, 1 )]-“, our analysis involves two oases: 0 < p < 1 and 1 < p < oo. 
Define 7 ( 5 , x) = as the lower inoomplete gamma funotion. The Laplace transforms 

of the distributions of and are given as follows. 

Lemma 3: The Laplaee transform of the distribution of aggregate reeeived power at the typieal 
SNq from the near PBs with gain Gm is 


where 


and 


H pM ( 5 ) — 
s,n ^ ' 




for 0 < p < 1 

(25a) 

GpM( 2 )(s), for 1 < p < cx). 

(25b) 

[p^ - exp {-sPpaGM)] | 

(26) 


GpMj2)(s) 



- exp {-sPpaGMp “) 


{sPpaG 


M 


7(1-^, sPpaGMj - 7 U - sPpaGMp~ 


(27) 


Proof: See Appendix A. 
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Lemma 4: The Laplace transform of the distribution of aggregate received power at the typical 
SNq from the far PBs with gain Gm is 

£pM(i)(s), for 0 < p < 1 (28a) 

£pM( 2 )(s), for 1 < p < cx), (28b) 

where 


£pM (i)(s) = exp ^ XpTirjf <j p^ - p^ exp (sPpaGM) - (sPpaGM) “ 7 ( 1 -^, sPpaGM 


and 


(29) 


'^p"( 2 )('S) =exp Xp7rr]f< exp [sPpaGMp ") 


- {sPpaGM)°‘l 1 -, sPpaGMp~ 

a 


(30) 


Proof: See Appendix B. ■ 

Based on Lemma 3 and Lemma 4, we obtain the Laplace transform of the distribution of Pg in 
the following proposition. 

Proposition 1: The Laplace transform of the distribution of aggregate received power at the 
typical SNq from all PBs under AD-WPT is given by 

C N N 

n (^) n ’ for 0 < p < 1 (31a) 


CpAs) = 


M=1 M=0 

N N 

n ^A^n( 2 ) (s) n M( 2 ) (s) , for 1 < p < CX). (31b) 

M=1 M=0 

As a special case of = 1, the Laplace transform of the distribution of aggregate received 
power at SNq from all PBs in omnidirectional WPT is given by 


Cpomni{s) = exp — Xp7i {sPpa) “ 7(1 -, sPpa 


(32) 


Proof: 


CpA-s) = E [exp {-sPs)] = E [exp (-s {Ps,n + -P^,/))] = E [exp {-sPs,n)] E [exp {-sPsj)] 


= E 


N 


exp -s PZ 


M=l 


E 


N 


exp [s^P, 


M 

s,f 


M=0 


N 


N 


n n (33) 


M=1 


M=0 


Substituting CpM (s) in Lemma 3 and CpM (s) in Lemma 4 into (|3^ . we obtain the Laplace 




transform given in Proposition 1. It is noted that CpA^) is continuous at p = 1. 
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C. Mean and Variance of Received Power 


In this subsection, we derive the closed-form mean and variance of the received power at SNq 
by taking the derivative of the Laplace transform in Proposition 1. These results will be useful 
in the approximation of the CCDF of received power in the next subsection and the average 
power maximization in Section IV. 

The average received power at SNq is given by 

E[P.] = —^ [log (>Cp.(s))] 1^=0- (34) 

This is the expectation of the aggregate received power at SNq from all PBs by taking over all 
possible location realizations of the PBs in spatial domain. By further derivations, the results 
are summarized in the following proposition. 

Proposition 2: At the typical SNq, the average received power in AD-WPT is given by 


^[Ps] = 


-/^p/\p(!T7r 


PpXpCTTT 


p2 (pp _ ^ 


[a 


1 — p 
- 


+ 


a 


{l-p^) ^ 2p 


2-a 


[a 


for 0 < p < 1 (35a) 


for 1 < p < oo, (35b) 


2)(l-p) ■ «-2 

where p is given in (fT^ and E(Ps) is continuous at p = 1. As a special case of = 1, the 
average received power at SNq in omnidirectional WPT is given by 

Ein” 




PpXpaira 


a — 2 


(36) 


Proof: See Appendix C. With N = 1, both (I35al) and (I35bl) equal (l3^ for all p. ■ 

In Proposition 2, for any given p, E[Ps] is increasing with the increased Pp, \p and N and 
is decreasing with the increased A^. Moreover, for any given set of {Pp, Ap, A^,, A^}, E,[Ps] is 
unimodal in p, i.e., there exists a unique p* that maximizes E[Ps], where E[Ps] is monotonically 
increasing for p < p* and is monotonically decreasing for p > p*. In Section IV, we will analyze 
the optimal p* that maximizes E[Ps]. 

The comparison of the average received power between AD-WPT and omnidirectional WPT 
is given in the following corollaries. 

Corollary 1: For 0 < p < oo, it follows that E[Ps] > E[P{”""*]. For p —)■ 0 and p —)■ cx), we 
have lS^[Ps] —)■ E[P°"*"'*]. 

Proof: See Appendix D. ■ 

From Corollary 1, we see that the average received power at SNq from all PBs in AD-WPT is 
higher than that in omnidirectional WPT. Next, we further discuss how the near PBs (||3fj|| < P) 
and far PBs (||3fi|| > p) influence the average received power. For comparison, we denote the 
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aggregate reeeived power from the PBs with ||Xj|| < p and ||Xj|| > p in omnidireetional WPT 
by and respeetively. 

Corollary 2: The ratio of the average reeeived power at SNq from the near PBs in AD-WPT 
and in omnidireetional WPT is given by 


E[P., 


■P 


N 


> 1 


(37) 


The ratio of the average reeeived power at SNq from the far PBs in AD-WPT and in omnidi¬ 
reetional WPT is given by 


mj] 

E[p;j™] 


(38) 


Proof: The proof is similar to that of Proposition 2 and thus omitted. ■ 

In Corollary 2, we see that the average reeeived power at SNq from the near PBs and far PBs in 
AD-WPT is greater than and equal to that in omnidireetional WPT as shown in (1371) and (l38l) . 
respeetively. The improvement of the average reeeived power from all PBs at SNq is thus due 
to the adaptive energy beamforming of the near PBs. 

By taking the seeond derivative of the Laplaee transform in Proposition 1, we obtain the 
varianee of the reeeived power at SNq, i.e., 


V[FJ = ^|log(£p.(s))]|,.„. (39) 

By further derivations, we summarize the results in the following proposition. 

Proposition 3: At the typieal SNq, the varianee of the reeeived power in AD-WPT is given 
by 


\pyT^< 


a 


P 


+ 


V[P«] = 


q; — 1 

A yV y-i 

^\m) [m-1 


y - 1 ) /+ 


a 


q; — 1 


X 


M =1 


for 0 < p < 1 (40a) 


I 




a — 1 


a — 1 

A yV y-I 
^ i mJ [m-1 


a — p‘^ p 2 2 q 


X 


M=1 


a — 1 


p^-MqM , for 1 < p < cx), (40b) 


where Y{Ps) is eontinuous at p = 1. As a speeial ease of = 1, the varianee of the reeeived 
power in omnidireetional WPT is 


v|p; 


XpPpa'^Tra 
a — 1 


(41) 
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Fig. 2. CCDF of the received power at SNq (Pp = 5 W, p = 2, Ap = 0.1, As = 0.2, i' = 0.1 m and a = 3). 


Proof: See Appendix E. With N = 1, both (I40al) and (I40bl) equal (ITO for all p. ■ 

In Proposition 3, V[P<j] is unimodal in p, i.e., V[Ps] is first inereasing and then deereasing 
with the inereased p. Given any p, V[Ps] is inereasing with Pp, Ap and N, and it is deereasing 
with As- We also eompare the varianee of the reeeived power at SNq between AD-WPT and 
omnidireetional WPT in the following eorollary. 

Corollary 3: For p —)■ 0 and p —>■ cxo, we have V[Ps] —)■ For 0 < p < cxo, it follows 

that V[P,] > V[p;"^™]. 

Proof: The proof is similar to that of Corollary 1 and thus omitted. ■ 

From Corollary 3, we see that varianee of reeeived power in AD-WPT is higher than that in 
omnidireetional WPT. Though AD-WPT improves the average reeeived power eompared with 
omnidireetional WPT, it also eauses more signifieant spatial power fluetuation. 

D. Characterization of Received Power via CCDF 

In this subseetion, we analyze the CCDF Fg of the reeeived power at the typieal SNq, whieh 
is the probability that Pg takes on a value greater than or equal to the threshold i.e., 

poo 

Fg = Pr (P, > Pf) = / f{Ps)dPg, (42) 

Jpth 

where f{Ps) is the probability density funetion (PDF) of Pg and ean be ealeulated from the 
inverse Faplaee transform of Cp^{s) in Proposition 1, i.e., 

f{Pg) = C-pl{s). 


(43) 
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In most cases, the direet derivation of the PDF from the inverse Laplaee transform is very 
ehallenging, if not possible, espeeially for the non-singular path loss model. In our problem, the 
new parameter p adds more diffieulty to the derivation. Even for omnidireetional transmission, 
the elosed-form PDF only exists for speeial ehoiees of parameters, e.g., a = 4 and singular path 
loss model 0. Studies have been shown that Gamma distribution gives a good fit to the power 
distribution for the homogeneous PPP lfT9l and heterogeneous PPP ll20l with non-singular path 
loss. In this work, we approximate the reeeived power distribution by Gamma distribution with 
seeond-order moment matching method, i.e., by matehing the mean and varianee of the two 
distributions, where the mean and varianee of Pg are given in Proposition 2 and Proposition 3, 
respeetively. Denote the Gamma funetion by r(fc) = The approximation of the 

CCDF of Pg is given in the following proposition. 

Proposition 4: At the typieal SNq, the approximated CCDF of the received power using 
Gamma distribution with seeond-order matehing is expressed as 


Fs 


1 - 


V{k) 


(44) 


where k = and 6 = are the shape parameter and seale parameter of the Gamma 

distribution, respeetively. 

Fig. 2 shows a good mateh between the simulation results of Fg and its approximation Fg. In 
some seenarios, an SN is aetive if the reeeived power is beyond the eonstant operational power 
threshold Pf". Then, the CCDF of Pg ean be regarded as the aetive probability of the SNs. It ean 
be proved that Fg in (l44l) is inereasing in E[Ps] and deereasing in V[P 5 ]. The inereased average 
received power and power fluetuation may improve or reduce the sensor aetive probability, 
respeetively. As shown in Proposition 2 and 3, both E,[Pg] and Y[Pg] first inerease and then 
deerease with the inereased p. We will further diseuss the above tradeoff and derive the optimal 
p* that maximizes the sensor aetive probability in Seetion V. For omnidireetional WPT, the 
CCDF of reeeived power and its approximation can be obtained by substituting 

(Ho]) into (021) and by substituting (l3^ and (l4TI) into (ITO . respeetively. 

Corollary 4: Fg inereases with the inereased Pp and/or Ap. 

Proof: It ean be proved that 6 and k are linear inereasing with Pp and Ap, respectively. As 
Pp inereases, 6 inereases and k remains a eonstant. Sinee Fg is an inereasing funetion of 6, it 
also inereases with the inereased of Pp. Similarly, as Ap inereases, k inereases and 6 remains a 
eonstant. Sinee Fg is an inereasing functions of k, it also increases with the inereased Ap. ■ 
Corollary 4 shows that inereasing the PB density or power improves the sensor aetive probability. 
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Fig. 3. Average received power versus charging radius (Pp = 10 W, i/ = 0.1 m, a = 3, Ap = 0.1 and N = 4). 


IV. Maximization of Average Received Power in AD-WPT 

In flexible-task WSN applieations, eaeh SN is assigned with flexible sensing tasks depending 
on the received energy, i.e., the SNs with high received power may handle more tasks than the 
low-power SNs for the benefit of the whole network. For example, the high-power SNs in a 
hierarchical network may work as cluster-heads |[I1 that collect data from the low-power SNs 
and coordinate sensing tasks among the SNs. The low-power SNs can also offload part of the 
computational processing tasks to the high-power SNs which have abundant resource If^ . To 
achieve energy efficient AD-WPT in flexible-task WSN, it is important to maximize the total 
received power over all SNs, which is equivalent to maximizing the average received power at 
the typical SNq. In this section, we design the optimal charging radius p* for maximizing E[Ps] 
in Proposition 2. 

In Fig. 3, we plot the average received power E[Ps] in AD-WPT, which outperforms the 
average received power in omnidirectional WPT for all p > 0. Moreover, E[Ps] first 

increases and then decreases with the increased p and there exists an optimal p* that maximizes 
E[Ps]. These results match well with Proposition 2 and Corollary 1 and are explained as follows. 
The received power at SNq from PBj depends on whether PBj radiates power towards SNq and 
the intensity of the radiated energy, which can be viewed as the power opportunity and power 
intensity of the PB, respectively. Both the power opportunity and power intensity are related 
to the number of beams of PBj and the distance ||Xj|| between PB* and SNq. We discuss the 
average received power at SNq from the near PBs (|| ViH < p) and far PBs (||XjH > p) as follows. 
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• Near PBs: The average reeeived power from the near PBs in AD-WPT is higher than that in 
omnidireetional WPT (see (iTTl) in Corollary 2). When the number of beams from a near PB 
deereases, the power intensity of this near PB inereases (see (|2])), and the power opportunity 
that the near PB radiates power towards SNq is with probability 1 sinee the near PB forms 
at least one beam towards SNq. 

. Far PBs: The average reeeived power from the far PBs in AD-WPT is the same as that 
in omnidireetional WPT (see (IMI) in Corollary 2). When the number of beams of a far 
PB deereases, the power intensity of this far inereases (see dU) and dH)), but the power 
opportunity that SNq reeeives energy from this far PBs deereases (see (II Sal) and (I18bl) l. 
and viee versa. From the average perspeetive, the effeets of power intensity and power 
opportunity of the far PBs eaneel with eaeh other. 

From the above diseussions, E,[Ps] outperforms mainly beeause of the high power 

intensity from the near PBs. We diseuss 'E[Ps] as follows. 

• As p —)■ 0, no SN is in the eharging regions. E[Ps] is equivalent to E[P°™™] sinee all PBs 
radiate power in N direetions. 

• As p inereases, a small number of SNs are ineluded in the eharging regions and the PBs that 
are elose to SNq beeome near PBs. When most PBs eoneentrate the transmit power from N 
beams into 1 beam, the power intensity is greatly enhaneed eompared with omnidireetional 
WPT. E[Ps] inereases with the inereased p due to the inereased number of the near PBs 
and inereased power intensity of the PBs. 

• As p further inereases, more seetors of the PBs are likely to be aetivated due to the inereased 
number of SNs in the eharging regions. When the number of beams of most PBs inereases 
from 1, 2, ■ • •, to N, the power intensity for eaeh beam deereases. There is a tradeoff 
between the further inereased number of the near PBs and the deereased power intensity. 
E[Ps] thus inereases and then deereases with the inereased p. 

• As p —)■ oo, all SNs are in the eharging regions and AD-WPT is again equivalent to 
omnidireetional WPT. 

In the following, we study the optimal eharging radius p* that maximizes E[Ps] in Proposition 
2 , i.e., 

PI ; E[Ps]* = max E[Ps]. (45) 

0<p<oo 

In Fig. 3, for omnidireetional WPT, E[P°”^"®] is regardless of A* whieh matehes with (l3^ . It 
is beeause eaeh PB radiates energy in all direetions without eatering to the loeations or density 
of the SNs. For AD-WPT, E[Ps] deereases with the inereased A^. With the inereased number of 
SNs in the eharging regions, PBs are more likely to radiate with more beams and less power 
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intensity, which thus reduces the average received power at SNq. In the following, we discuss 
the optimal p* for different under AD-WPT. 

Case 1: Low SN Density Regime. When the SN network density is low, e.g., = 0.2 in Fig. 

3, we have dp ^ lp=l ^ Since (I35bl) is unimodal in p and (I35al) is an increasing function 
of p, the optimal charging radius p* G (1, cxo) is the stationary point of (I35bl) . Taking the first 
derivative of (I35bl) with respect to p, we have 


a[E[P,]|i<,<oo] 


dp 


2 iPp Ap TT (T 


pi « p{a-2p^ °) 

l-p (a-2 )(l-p) 


XgT^pp (a — 2p^ “) (l — 
(1 -p)2(a -2)iV 


(46) 


The optimal charging radius p* is the unique solution of = 0. Though p* is not 

analytically tractable, we can search it numerically using one-dimensional searching method. 
Case 2: Medium SN Density Regime. When the SN network has a medium density, e.g., A^ = 0.8 
in Fig. 3, we have = 0. In this case, (I35al) is an increasing function of p and (I35bl) is 

a decreasing function of p. The optimal charging region radius is at the point of p* = 1. 

Case 3: High SN Density Regime. When the SN network has a high density, e.g., A^ = 1.6 in 
Fig. 3, we have < 0. In this case, (I35al) is unimodal in p and (I35bl) is a decreasing 

function of p. The optimal charging radius p* G (0,1) is the stationary point of (I35al) . Taking 
the first derivative of (I35al) . we have 


a[E[P,]|o<,<i] 

dp 


2PpXp7ra 


XsTTp^pip^ -p) p{p-p^) + Xs'Kp^p {p^ ^ 

N{l-pf ^ l-p 



(47) 


The optimal charging radius p* is the unique solution to = 0. Similar to Case 1, p* 

is not analytically tractable but can be searched numerically. 

It can be proved that and are of the same sign at the point of p = 1. 

The procedure of obtaining the optimal p* is summarized in Algorithm 1. More numerical results 
will be shown in Section VI-A. 


V. Maximization of Sensor Active Probability in AD-WPT 

In the previous section, we discussed the optimal AD-WPT design in flexible-task WSN 
scenario where the energy consumption levels or tasks vary for different SNs. In some other 
scenarios, e.g., environmental measurement ll22l and surveillance monitoring [|23l systems, the 
sensing information from each SN is equally important and mutually irreplaceable. For example, 
in a forest fire detection systems |[I1, SNs are randomly deployed in a forest collecting temperature 
and humidity data independently. In these scenarios, the SNs are assigned with equal sensing 
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Algorithm 1 Solving the optimal charging radius in PI: 

1: 

Calculate Di(p) - 

and D 2 {p) 

_ 5[lE[Ps]|l<p<oo] 

dp 

2 : 

if either Di{p = 1) < 0 or D 2 

A 

0 

then 

3: 

p* is the solution to Di{p) = 

= 0 



4: 

else if either Di {p = 1) = 0 or D 2 {p = 

1) 

= 0 then 

5: 

p* = l 




6: 

else if either Di {p = 1) > 0 or D 2 {p = 

1) 

> 0 then 

7: 

p* is the solution to T> 2 (p) = 

= 0 



8 : 

end if 





tasks with a minimum operational power requirement [[24l. i.e., an SN is active if the received 
power is beyond the target energy threshold. To achieve higher sensing diversity, it is important to 
allow more SNs operating with sufficient power. In this section, we analyze the optimal charging 
radius p* in AD-WPT to maximize the active probability Fg of the SNs. 

As discussed in Section IV, the decreased number of beams at the PBs improves the radiated 
power intensity, which enhances the average received power at SNq in AD-WPT compared with 
omnidirectional WPT. However, the decreased number of beams may not enhance the sensor 
active probability Fg due to the interplay between the power intensity and power opportunity. 

• Near PBs: The near PBs help improve the sensor active probability in AD-WPT compared 
with that in omnidirectional WPT. Since the near PBs always radiate energy towards SNq 
with probability 1 (see (fTSl) ! and antenna gain greater than 1 (see dH)), the received power 
from the near PBs in AD-WPT is higher than that in omnidirectional WPT. With the de¬ 
creased number of beams from the near PBs, the power intensity increases, which increases 
the received power from the near PBs and may improve the sensor active probability in 
AD-WPT. 

• Far PBs: The far PBs can reduce the sensor active probability in AD-WPT compared with 
that in omnidirectional WPT. Since the far PBs may not radiate energy towards SNq, the 
received power from a far PB in AD-WPT is higher than that in omnidirectional WPT or 
zero if SNq is inside or outside the beamforming directions of the PB, respectively. With the 
decreased number of beams from the far PBs, the power intensity of the far PBs increases 
(see (HI) and dH)), but the power opportunity to receive energy from the far PBs at SNq 
decreases (see (IISal) and (I18bl) l. Since SNq has a higher chance to fall outside the radiation 
directions of the far PBs, the received power from the far PBs is more likely to decrease, 
which may reduce the sensor active probability in AD-WPT. 
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Fig. 4. Sensor active probability versus charging radius with various PB power (Ap = 0.1, As = 0.2, Pl^ = 0.1 mW, u = 0.1 
m and a = 3). 


As p increases from 0 to oo, the number of beams of most PBs deereases from N to 1, and 
then inereases from 1, 2, • • •, to With fewer beams, the inereased power intensity of the 
near PBs and deereased power opportunity of the far PBs have positive and negative impaets on 
the sensor aetive probability Fg, respeetively. Whether the near PBs or the far PBs dominate Fg 
depends on the PB power/density and radius p. If the PB power/density is low or p is large, the 
near PBs dominate Fg due to the severe power attenuation of the far PBs, and viee versa. 

In the following, we analyze the optimal eharging radius p* that maximizes the sensor aetive 
probability Fg at the typieal SNq, i.e., 

P2 : F* = max Fg. (48) 

0<p<oo 

The simulation results of the sensor aetive probabilities Fg in AD-WPT and in omnidi- 

reetional WPT are plotted against p with various Pp in Fig. 4. As shown in Corollary 4, the 
sensor aetive probability inereases with inereased PB power Pp or PB density Xp. We diseuss 
Fg by eonsidering different power regimes of the PBs. 

Case 1: Low PB power/density regime. When the PBs have low power Pp and/or density Xp, 
e.g., Pp = 1 W in Fig. 4, the near PBs dominate Fg and we have Fg > Fg inereases and 

then deereases with p mainly due to the deereased and inereased power intensity of the PBs, 
respeetively. There exists an optimal p* (e.g., p = 1.5) that maximizes Fg. 

Case 2: Medium PB power/density regime. When the PBs have medium power Pp and/or density 
Xp, e.g., Pp = 3 W in Fig. 4, the far PBs dominate Fg for small p with Fg < F^™, and the 
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near PBs dominate Fg for large p with Fg > respectively. There exists an optimal p* 

(e.g., p = 2.25) that maximizes Fg in the region of Fg > 

Case 3: High PB power/density regime. When the PBs have high power Pp and/or density 
Ap, e.g., Pp = 10 W in Fig. 4, we have Fg < Due to the high PB power, the sensor 

active probability for omnidirectional WPT is high. For AD-WPT, Fg decreases and then 

increases with p mainly due to the decreased and increased power opportunity of the far PBs, 
respectively. As p —)■ 0 or p ^ cxd, we have Fg —)■ F°'^'^\ The optimal p* is approaching 0 or 
oo as in omnidirectional WPT. 

From the above discussions, we see that the maximized sensor active probability in AD-WPT 
with the proper selection of charging radius p* is larger than or at least equivalent to that in 
omnidirectional WPT. As discussed in Section III-D, Fg is not analytically tractable but can be 
well approximated using Gamma distribution. In the following, we analyze the optimal charging 
radius p* that maximizes the approximation of sensor active probability Fg in Proposition 4, i.e.. 


P3 ; F* = max Fg. 

0<p<oo 


(49) 


To solve the one-dimensional problem in P3, the optimal radius is one of the stationary points 
of Fg. Taking the derivative of Fg with respect to p yields 


dp 




V ^ ) 


™V[P.]+E[P,]^ 


dp 


|v[Fj]^r(fc) 


a[E[P. 

dp 


-V[P,] + [E[P,]] 


2 av[p.] 

dp 


p _ § _ 

/o‘ * 




\n{t)dt 


|v[P.ll"r(i) 


I 1 pF 

7 


a[E[Pa: 




IVIPJI" |r (A^)]" 


(50) 


where dW[Pg]/dp is given in Section IV and dY[Ps]/dp can be obtained via similar approaches. 
The number of stationary points of Fg is less than or equal to two since dFg/dp = 0 has at most 
two solutions for p G (0, cx)). The procedure to obtain the optimal radius for P3 is summarized 
in Algorithm 2. 


VI. Numerical Results 

In this section, we present the simulation results of the maximized average received power 
for flexible-task WSN and the maximized sensor active probability for equal-task WSN under 
the proposed AD-WPT scheme, respectively. The performance of omnidirectional WPT scheme 
is used as a comparison benchmark. Throughout this section, we set a = —41.9842 dB, where 
the wavelength is z/ = 0.1 m and reference distance is do = 1 
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Algorithm 2 Solving the optimal charging radius in P3: 

1 : Find J as the number of stationary points of F* for p G (0, cx)) and calculate ^|p=o+ 

2 : if J = 1 and ^|p=o+ < 0 (as in Case 1) then 

3: p* is the single stationary point of Fs 

4 : else if J = 2 and ^|p=o+ > 0 (as in Case 2) then 

5: p* is the stationary point of Fg with larger value of p 

6: else if J = 1 and ^|p=o+ > 0 (as in Case 3) then 

7: p* is approaching 0 or cxo 

8 : end if 



(a) (b) 


Fig. 5. (a) Optimal charging radius p* for average power maximization versus N (Ap = 0.1, As = 0.2 and a = 3). (b) 

Maximized average received power E[Ps]* versus N (Ap = 0.1, As = 0.2 and a — 3). 


A. Maximized Average Received Power for Flexible-task WSN 

Fig. 5 (a) and Fig. 5 (b) show that both the maximized average received power E[Ps]* and the 
corresponding optimal charging radius p* increase with the increased number of PB sectors N. 
As N increases, the PBs are able to form narrower energy beams with higher power intensity 
towards the intended SNs. As a result, the coverage of PBs in AD-WPT extends and it is 
more beneficial to use a larger charging radius p* as shown in Fig. 5 (a) to serve more SNs 
efficiently. With the decreased beamwidth, the power intensity of the near PBs increases, which 
thus improves E[Ps]* as shown in Fig. 5 (b). This is similar to the effect of decreasing the 
number of beams of the PBs as discussed in Section IV. 
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As, SN density 
(a) 


As, SN density 

(b) 


Fig. 6. (a) Optimal charging radius p* for average power maximization versus As (Ap = 0.1, a = 3 and N = 4). (b) 

Maximized average received power E[Ps]* versus As (Ap = 0.1, a = 3 and N = 4). 


Fig. 6 (a) and Fig. 6 (b) illustrate that the maximized average received power E[Ps]* and 
the corresponding optimal charging radius p* decrease with the increased SN density A^. As 
As increases, more sectors of PBs are activated. The PBs form more energy beams with lower 
power intensity towards the intended SNs. As a result, the PBs shrink the charging radius in 
AD-WPT as shown in Fig. 6 (a) to serve fewer SNs efficiently. In Fig. 6 (b), E[Ps]* decreases 
with the increased A^ due to the decreased power intensity of the near PBs. 

From Fig. 5 (a) to Fig. 6 (b), we observe that E[Ps]* increases linearly with the increased PB 
power Pp, but p* is regardless of Pp. It can also be deduced that increasing Ap has a similar 
impact on E[Ps]* as increasing Pp. Moreover, when Pp is relatively high (e.g., Pp = 8 W), 
increasing N or decreasing A^ causes more significant improvement of E[Ps]* as shown in Fig. 
5 (b) and Fig. 6 (b), respectively. 

B. Maximized Sensor Active Probability for Equal-task WSN 

In Fig. 7 (a), the maximized sensor active probability F* increases with the increased N. As 
discussed in section V, the improvement of the sensor active probability in AD-WPT compared 
with omnidirectional WPT is mainly due to the high power intensity of the near PBs. As N 
increases, the beamwidth of the PBs decreases, which improves the power intensity and thus 
improves F*. Fig. 7 (b) shows the maximized sensor active probability F* decreases with the 
increased A^. As A^ increases, the near PBs are likely to form more beams with lower power 
intensity, which therefore reduces P*. Furthermore, we notice that increasing N or decreasing 
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Fig. 7. (a) Maximized sensor active probability F* versus N (Ap = 0.1, As = 0.2, = 0.1 mW and a = 3). (b) Maximized 

sensor active probability F* versus As (Ap = 0.1, A'’ = 4, P*^ = 0.1 mW and a = 3). 


Xg causes more signifieant improvement of F* for relative small Pp (e.g., Pp = 2 W). The 
improvement is less signifieant for relatively large PB power (e.g., Pp = 8 W) sinee the sensor 
aetive probability in omnidireetional WPT is already high and may not be mueh improved by 
AD-WPT. 

C. Comparison with Other Power Allocation Schemes 

For the AD-WPT seheme in Seetion II-A, we adopt uniform power alloeation for the PBs, 
i.e., uniformly alloeating the PB power among all aetive seetors that have at least one SN. If 
the exaet number of SNs in eaeh seetor is known, the PBs may adopt unequal power alloeation 
sehemes whieh alloeate the PB power aeeording to the number of SNs in eaeh seetor. In this 
subseetion, we mainly diseuss two other power alloeation sehemes: greedy seheme and robust 
seheme. In greedy seheme, eaeh PB alloeates all power to the seetor that has the largest number 
of SNs and no power to all other seetors. It ean be easily shown that this seheme provides the 
maximum sum reeeived power of all SNs in the eharging region of a PB. In robust seheme, 
eaeh PB alloeates power proportionally to the number of SNs in eaeh seetor. 

In Fig. 8 (a) and Fig. 8 (b), we eompare the maximized average reeeived power and the 
maximized sensor aetive probability for the three power alloeation sehemes, respeetively. Fig. 8 
(a) shows that the maximized average reeeived power for the three sehemes is similar with very 
minor gaps. Greedy seheme performs the best and robust seheme slightly outperforms uniform 
seheme in terms of average reeeived power. In Fig. 8 (b), we see that the maximized sensor 















25 



Fig. 8. (a) Comparison of the maximized average received power for the three power allocation schemes (Ap = 0.1, As = 0.2, 

Af = 4 and a = 3). (b) Comparison of maximized sensor active probability for the three power allocation schemes (Ap = 0.1, 
As = 0.2, Pi'* = 0.1 mW, AT = 4 and a = 3). 


active probability of robust scheme is the highest, which slightly outperforms that of uniform 
scheme. The sensor active probability of greedy scheme is the lowest. 

To sum up, greedy scheme shows the highest average received power but at the cost of lowest 
sensor active probability. It is because the single beam strategy in greedy scheme improves the 
power intensity but also reduces the power opportunity towards SNq. Moreover, robust scheme 
outperforms uniform scheme in both average power and active probability of the SNs, but the 
improvement is insignificant. For greedy and robust schemes, the exact number of SNs in each 
sector is required and the derivation of the distribution of the received power in a heterogeneous 
network is more complicated than uniform scheme since the gain of the PB becomes a continuous 
instead of discrete variable. For uniform scheme, each PB needs only the information of the 
existence of SNs in each sector. It provides acceptable average power and active probability 
with less implementation complexity. 


VII. Conclusions 

In this paper, we proposed an AD-WPT scheme in a large-scale sensor network, where the PBs 
charge the SNs by adapting the energy beamforming strategies to the nearby SN locations. By 
using stochastic geometry, we derived the closed-form expressions of the distribution metrics, 
i.e., Laplace transform, mean and variance, of the aggregate received power at a typical SN. 
The approximation of the CCDF of received power is obtained using Gamma distribution with 
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second-order moment matching. For flexible-task and equal-task WSN, the optimal radii in 
AD-WPT were designed to maximize the average received power and maximize the sensor 
active probability, respectively. The results show that the maximized average received power and 
maximized sensor active probability increase with the increased density and power of the PBs, 
while they decrease with the increased density of the SNs and energy beamwidth. Moreover, the 
optimal AD-WPT is more energy efficient than omnidirectional WPT by achieving equivalent 
average received power or sensor active probability with less transmit power consumptions. 


Appendix 

A. Proof of Lemma 3 

Taking the Laplace transform of (|2^ . we have 




exp 


-sPpCr ^ Gm [max(||Xi|l, !)]■ 




= E 


= exp 


exp [-sPpaGu [max (||Xi|l, 1)] "] 

27r np 




I I [l - exp [-sPpCrGju [max (r, 1)] rdf dr . (51) 

Jo Jo 

The last step is obtained by applying probability generating functional (PGFL) Proposition 
2.12], where r and ip denote the radial coordinate and angular coordinate in polar coordinate 
system. For 0 < p < 1, (|5TI) is further derived as 


(^) = exp 


27iXpT]^ I [1 — exp (—sPpCtGm)] rdr 


(52) 


For 1 < p < oo, (|5TI) is further derived as 


= exp 


—27rApp, 


M 


Uo 


[1 — exp (—sPp(jGm)] pdr 


[l — exp (y—sPpaGMP “)] rdr 


(53) 


From (15^ and (l5^ . we can easily obtain (|2^ and (l27l) in Lemma 3, respectively. 


B. Proof of Lemma 4 

Taking the Laplace transform of (|2^ . we have 

PpM (s) = E [exp (-sP,j)] 


^ 27 r poo 


= exp 


Kvf 


M 


[l — exp [—sPpCrGM [max (r, 1)] “]] rdf dr 
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(54) 
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For 0 < p < 1, (l54l) is derived as 


£pM (s) = exp— 27rApPj j [1 — exp {—sPpaGM)]rdr 


— 27rApp|^ J [l — exp [—sPpaGMr~'^)] rdr ^ 
For 1 < p < cxD, (l54l) is derived as 

CpM (s) = exp < — 27rApP^^ f [l — exp (^—sPpaGMr~°‘)] rdr 


(55) 


(56) 


From (1551) and (1561) . we ean easily obtain (1291) and (1301) in Lemma 4, respeetively. 

C. Proof of Proposition 2 

In this appendix, we derive E(Ps) in Proposition 2. 

i) 0 < p < 1; Taking the first derivative of the Laplaee transform in Proposition 1 for 
0 < p < 1, we have 

E[Pj]|o<p<l = [log ('(^-Ps('S)|o<p<l)] \s=0 


^ d ^ d 

-ApTTp^ ^ — [exp {-sPpaGM)] |s=o + ApVrp^ ^ ds {-sPp(tGm)] |s=o 
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+ XpTT Vj 

M =0 

By further derivation, we have 

d 

ds 


7 ( 1 - 


U=o- 
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M 


and 


Substituting (1581) and (1591) into (1571) yields 
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We further obtain 
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Substituting (IMI) and (l62l) baek into (l60l) . we obtain (135al) in Proposition 2. 
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2j 1 < p < oo; Taking the first derivative of the Laplaee transform in Proposition 1 for 
1 < p < oo, we have 

E[P.]|i<p<oo = -^[log(£ r’s('5) |i<p<oo)] |s=o 

' N N 

vfGu - ^ 

M=0 M=1 

Substituting (IMl) and (1^ into (1^ . we obtain (I35bl) in Proposition 2. 


PpXp7T(T 
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D. Proof of Corollary 1 

We compare (l3^ with (I35al) and (I35bl) . respectively. Firstly, (13Sal) is compared with (l3^ . 

_ ^ 

E [Ps] |o<p<i - E [Pr"*] = ^pApVrap^^^, (64) 

where p was given in (fT^ . For p —)■ 0, it has p —)■ 1 and E [P^] |o<p<i —)■ E[P°™™]. For 
0 < p < 1, it has < p < 1 which leads to E [P^] |o<p<i > E [P°”^”*]. 

Secondly, (I35bl) is compared with (l3^ . 

E [P.] |i<p<oo - E (Pr"‘) = PpAp7ra^^^4-T^- 

a — 2 1—p 

For p —)■ oo, it has p —)■ 0 and E [P^] |i<p<oo —^ E For 1 < p < oo, it has p'^ < p < 1. 

Since a > 2, we further have p^“" < 1 and a > 2p^““. Then, it is proved that E [P*] |i<p<oo > 
E[P°mni^_ 


E. Proof of Proposition 3 

Taking the second derivative of the Laplace transform in Proposition 1, we have 
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By further derivation, we have 
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Substituting (1671) and (l6^ into (1^ yields (I40al) . Similarly, we can obtain V[Ps]|i<p<oo in (I40bl) . 
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